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With the aid of the quasiclassical Eilenberger formalism we develop a theory of non-local electron 
transport across three-terminal ballistic normal-superconducting-normal (NSN) devices with spin- 
active NS interfaces. The phenomenon of crossed Andreev reflection (CAR) is known to play the key 
role in such transport. We demonstrate that CAR is highly sensitive to electron spins and yields a 
rich variety of properties of non-local conductance which we describe non-perturbatively at arbitrary 
voltages, temperature, spin-dependent interface transmissions and their polarizations. Our results 
can be applied to multi-terminal hybrid structures with normal, ferromagnetic and half-metallic 
electrodes and can be directly tested in future experiments. 
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I. INTRODUCTION 

Low energy electron transport in hybrid structures 
composed of a normal metal (N) and a superconduc- 
tor (S) is governed by Andreev reflection 1 . (AR) which 
causes non-zero subgap conductance^ of such structures. 
AR remains essentially a local effect provided there ex- 
ists only one NS interface in the system or, else, if the 
distance between different NS interfaces greatly exceeds 
the superconducting coherence length £. If, however, 
the distance L between two adjacent NS interfaces (i.e. 
the superconductor size) is smaller than (or compara- 
ble with) £, two additional non-local processes come into 
play (see Fig. 1). One such process corresponds to direct 
electron transfer between two N-metals through a super- 
conductor. Another process is the so-called crossed An- 
dreev reflection 3 -" 4 . (CAR): An electron penetrating into 
the superconductor from the first N-terminal may form 
a Cooper pair together with another electron from the 
second N-terminal. In this case a hole will go into the 
second N- metal and AR becomes a non-local effect. Both 
these processes contribute to the non-local conductance 
of hybrid multi-terminal structures which has been di- 
rectly measured in several recent experiments 5 -!^. 

Theoretically the non-local conductance of NSN hy- 
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FIG. 1: Two elementary processes contributing to non-local 
conductance of an NSN device: (1) direct electron transfer 
and (2) crossed Andreev reflection. 



brids was analyzed within the perturbation theory in 
the transmission of NS interfaces in Refs. [H0 where it 
was demonstrated that in the lowest order in the inter- 
face transmission and at T = CAR contribution to 
cross-terminal conductance is exactly canceled by that 
from elastic electron cotunneling (EC), i.e. the non-local 
conductance vanishes in this limit. Thus, in order to 
determine the scale of the effect it is necessary to in- 
clude higher order (in the barrier transmission) terms 
into consideration. The corresponding analysis was em- 
ployed in Refs. Holllll by means of effective "dressing" 
of both EC and CAR contributions by higher order pro- 
cesses and more recently by the present authors^ with 
the framework of quasiclassical formalism of Eilenberger 
equations. We note that the results 1 ^ and 1 ^ disagree be- 
yond perturbation theory since "dressing" procedure 10 
does not account for all higher order processes. Hence, 
the approach^ is in general insufficient to correctly de- 
scribe non-trivial interplay between normal reflection, 
tunneling, local AR and CAR to all orders in the inter- 
face transmissions. We will return to this issue further 
below. 

Another interesting issue is the effect of disorder. It 
is well known that disorder enhances interference effects 
and, hence, can strongly modify local subgap conduc- 
tance of NS interfaces in the low energy limi t 13 : 14 ' 15 . 
Non-local conductance of multi-terminal hybrid NSN 
structures in the presence of disorder was recently stud- 
ied in Refs. [IallWMl9l . Brinkman and Golubov made 
use of the quasiclassical formalism of Usadel equations 
and proceeded perturbatively in the interface transmis- 
sions. Duhot and Melin 18 discussed the impact of weak- 
localization-type of effects inside the superconductor on 
non-local electron transport in NSN structures. Morten 
et a/Ji employed the circuit theory (thereby going be- 
yond perturbation theory in tunneling) and considered 
a device with normal terminals attached to a supercon- 
ductor via an additional normal island (dot)^ -. Very 
recently a similar structure with a superconducting dot 
was analyzed 1 ^ providing a rather general theoretical 
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framework to study non-local electron transport in multi- 
terminal NSN structures in the presence of disorder and 
non-equilibrium effects. 

Yet another interesting subject is an interplay between 
CAR and Coulomb interaction. The effect of electron- 
electron interactions on AR was investigated in a number 
of paper s 15 ' 21 ' 22 . Interactions should also affect CAR, 
e.g., by lifting the exact cancellation of EC and CAR 
contributions 2 ^ already in the lowest order in tunneling. 
A complete theory of non-local transport in realistic NSN 
systems should include both disorder and interactions 
which remains an important task for future investiga- 
tions. 

An important property of both AR and CAR is that 
these processes should be sensitive to magnetic properties 
of normal electrodes because these processes essentially 
depend on spins of scattered electrons. One possible 
way to demonstrate spin-resolved CAR is to use ferro- 
magnets (F) instead of normal electrodes ferromagnet- 
superconductor-ferromagnet (FSF) structure o 5 ' 10 ' 24 ' 25 . 
First experiments on such FSF structures^ illustrated 
this point by demonstrating the dependence of non-local 
conductance on the polarization of ferromagnetic termi- 
nals. Hence, for better understanding of non-local effects 
in multi-terminal hybrid proximity structures it is de- 
sirable to construct a theory of spin-resolved non-local 
transport. In the lowest order in tunneling this task 
was accomplished in Ref. [|| For FSF structures higher 
orders in the interface transmissions were considered in 
Refs^i. 

In this paper we are going to generalize our quasi- 
classical approach-^ and construct a theory of non-local 
electron transport in ballistic NSN structures with spin- 
active interfaces to all orders in their transmissions. Our 
model allows to distinguish spin-dependent contributions 
to the non-local conductance and to effectively mimic 
the situation of ferromagnetic and/or half- metallic elec- 
trodes. 

The structure of the paper is as follows. In Sec. 2 
we introduce our model and discuss the quasiclassical 
formalism supplemented by the boundary conditions for 
Green-Keldysh functions which account for electron scat- 
tering at spin-active interfaces. In Sec. 3 we employ this 
formalism and develop a theory of non-local spin-resolved 
electron transport in NSN structures with spin-active in- 
terfaces. Our main conclusions are briefly summarized in 
Sec. 4. 



II. THE MODEL AND FORMALISM 

Let us consider three-terminal NSN structure depicted 
in Fig. [5J We will assume that all three metallic elec- 
trodes are non-magnetic and ballistic, i.e. the elec- 
tron elastic mean free path in each metal is larger than 
any other relevant size scale. In order to resolve spin- 
dependent effects we will assume that both NS interfaces 
are spin-active, i.e. we will distinguish "spin-up" and 




FIG. 2: Schematics of our NSN device. 



"spin-down" transmissions of the first (Dij and -Dix) and 
the second (-D 2 | an d D 2 i) SN interface. All these four 
transmissions may take any value from zero to one. We 
also introduce the angle tp between polarizations of two 
interfaces which can take any value between and 2tt. 

In what follows effective cross-sections of the two in- 
terfaces will be denoted respectively as Ai and A 2 - The 
distance between these interfaces L as well as other ge- 
ometric parameters are assumed to be much larger than 
y/Ai,2, i- e - effectively both contacts are metallic con- 
strictions. In this case the voltage drops only across SN 
interfaces and not inside large metallic electrodes. Hence, 
nonequilibrium (e.g. charge imbalance) effects related to 
the electric field penetration into the S-electrode can be 
neglected. In our analysis we will also disregard Coulomb 
effects! 5 - 2 ^. 

For convenience, we will set the electric potential of 
the S-electrode equal to zero, V = 0. In the presence of 
bias voltages V\ and V 2 applied to two normal electrodes 
(see Fig. ^ the currents I\ and 1% will flow through SNi 
and SN2 interfaces. These currents can be evaluated with 
the aid of the quasiclassical formalism of nonequilibrium 
Green-Eilenberger-Keldysh functions 2 ^ g R < A > K which we 
briefly specify below. 



A. Quasiclassical equations 

In the ballistic limit the corresponding Eilenberger 
equations take the form 



eh + eV(r, t) - A(r, t), g R ' A ' K (p F ,e, r, t) 



+iv F Vg R > A > K (p F ,e,r,t) = 0, 



(1) 



where [a, b] = ab— ba, e is the quasiparticle energy, p F = 
mv F is the electron Fermi momentum vector and T3 is the 
Pauli matrix in Nambu space. The functions g R ' A ' K also 
obey the normalization conditions (g R ) 2 = {g A ) 2 = 1 
and g R g K + g K g A = 0. Here and below the product of 
matrices is defined as time convolution. 

Green functions g R ' A > K and A are 4x4 matrices in 
Nambu and spin spaces. In Nambu space they can be 
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parameterized as 



, n R.A,K fR,A,K 
' — I fR,A,K gR.A.K ) » ~ 



Ai<7 2 

A*i<7 2 



( 2 ) 

where g R ' A ' K , jR,a,k ^ Jr,a,k ^ gR,A,K are 2 x 2 matrices 
in the spin space, A is the BCS order parameter and <7i 
are Pauli matrices. For simplicity we will only consider 
the case of spin-singlet isotropic pairing in the supercon- 
ducting electrode. The current density is related to the 
Keldysh function g K according to the standard relation 

j( r >*) = --jp J de(v F Sp[T 3 g K (p F ,e,r,t)}) , (3) 

where N = mp F /2ir 2 is the density of state at the Fermi 
level and angular brackets (...) denote averaging over the 
Fermi momentum. 



B. Riccati parameterization 

The above matrix Green-Keldysh functions can be con- 
veniently parameterized by four Riccati amplitudes <y R ' A : 
^R,A an( j £ WQ "distribution functions" x K , x K (here and 
below we chose to follow the notations 30 ): 



• A, 
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FIG. 3: Riccati amplitudes for incoming and outgoing trajec- 
tories from the both sides of the interface. 



Depending on the particular trajectory it is also conve- 
nient to introduce a "replica" of both Riccati amplitudes 
and distribution functions which - again following the 
notation s 30 ! 31 - will be denoted by capital letters T and 
X. These "capital" Riccati amplitudes and distribution 
functions obey the same equations (f7|)- (fl0|) with the re- 
placement 7 — > r and x — ► X. The distinction between 
different Riccati amplitudes and distribution functions 
will be made explicit below. 



xK _ ~ R / x K - ~f R x K ^ A - 7 R x K + X K 1 A \ ~ A 



-^ R X K + X K J A X K — J R X K J A 



m/1 + 7 M 7 M 2 7 ^ , 

_2^ R ,A _i _ ^R,A^R,A I \°) 



where functions j R > A and j R ' A are Riccati amplitudes 

g R > A = 4 .V 

and N R ' A are the following matrices 
N R > A 



(4) 



{l _ 1 R,Az j R,A)-l . 

(1 _ !yR,A 7 R,Ayl I ■ W 







With the aid of the above parameterization one can iden- 
tically transform the quasiclassical equations fl} into the 
following set of effectively decoupled equations for Ric- 
cati amplitudes and distribution functions^ - 



iv F V^ R ' A + [e + eV{r, t)]j R ' A + j R ' A [e - eV(r, t)} 

= j R > A A*ia 2 j R > A -Aia 2 , 

iv F Vj R > A - [e - eV(r, t)]^ R ' A - j R ' A [e + eV(r, t)] 

= j R > A Aia 2 j R > A - A*ia 2 , 



(J) 



iv F Vx K +[e + eV{r, t)]x K - x K [e + eV(r, t)} 
~j R A*ia 2 x K - x K Aia 2 j A = 0, 



iv F \7x K - [e - eV(r, t)]x K + x K [e - eV(r, t)) 



z.K , ~K\ 



^ R Aia 2 x K - i K A*ia 2l A = 0. 



Z.K ~K 



(8) 
(9) 



(10) 



C. Boundary conditions 

Quasiclassical equations should be supplemented by 
appropriate boundary conditions at metallic interfaces. 
In the case of specularly reflecting spin-degenerate inter- 
faces these conditions were derived by Zaitsev^i and later 
generalized to spin-active interfaces in Ref. I32I . 

Before specifying these conditions it is important to 
emphasize that the applicability of the Eilenberger quasi- 
classical formalism with appropriate boundary conditions 
to hybrid structures with two or more barriers is, in gen- 
eral, a non-trivial issu o 28 i 29 . Electrons scattered at differ- 
ent barriers interfere and form bound states (resonances) 
which cannot be correctly described within the quasiclas- 
sical formalism employing Zaitsev boundary conditions 
or their direct generalization. Here we avoid this prob- 
lem by choosing the appropriate geometry of our NSN 
device, see Fig. 2. In our system any relevant trajectory 
reaches each NS interface only once whereas the proba- 
bility of multiple reflections at both interfaces is small in 
the parameter A\A 2 /L A <C 1. Hence, resonances formed 
by multiply reflected electron waves can be neglected, 
and our formalism remains adequate for the problem in 
question. 

It will be convenient for us to formulate the boundary 
conditions directly in terms of Riccati amplitudes and 
the distribution functions. Let us consider the first NS 
interface and explicitly specify the relations between Ric- 
cati amplitudes and distribution functions for incoming 
and outgoing trajectories, see Fig. [3] The boundary con- 
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ditions for rf , ly and X-f can be written in the form 31 

If = rfo? S++ tf^Sf,,, (11) 

r? = S u tfr? r + Soffit, (12) 

X\ = Tu^^ir + tli x Tjtlr ~ a ll^V ^lr- (13) 

Here we denned the transmission (i), reflection (r), and 
branch-conversion (a) amplitudes as: 

r£ = +[(^i)"^i + i - fer^irWi)- 1 , (14) 
*S = - (^lO -1 ^]- 1 ^!')- 1 . (15) 

ffr = - Sii'^iO" 1 ] -1 , (16) 

& = -O^l')- 1 ^!!^!)- 1 ~ Sn^!,)- 1 ]- 1 , (17) 

(18) 
(19) 



4 = (rf i 5 11 -5 11 7?)(/9f 1 0- 1 , 



where 



ffi = St 3 ~ l?S±lf; 0$ = Sji - ifs^?, (20) 



m = Sij - ifS^f, 13% = S± - itSffif- (21) 

Similarly, the boundary conditions for ff , rf , and 
take the form: 

rf = f^S ll +if l ^S vx , (22) 

rt = St 1 rfr? r + S+ l7 f,tt (23) 

X? = r&frl + i%x$t£. - afaZat (24) 

where 

r?t = +[(^i)- 1 S 11 - fer^J-^i)" 1 , (25) 

*S = -[(^J-^u - (ffiv)- 1 ^}- 1 ^,)- 1 , (26) 

= +(^i)" 1 |^iC^i)- 1 - SUPvv)- 1 }- 1 , (27) 

& - HttvY^SiMi)- 1 - S+Mv)' 1 ]- 1 , (28) 

af/ = (rf5+-5+7f)(/3f 1 0- 1 , (29) 

(30) 



afr = (P?y)- 1 (Si 1 r?-rfS 11 ) 



Boundary conditions for Ty' A , Fj)'^ 1 , X" and Xy can be 
obtained from the above equations simply by replacing 
1<->1'. 

The matrices 5n, Sii', Si/i, and 5i'i' constitute the 
components of the 5-matrix describing electron scatter- 
ing at the first interface: 



S = 



S\i S\y 

Sin Syy 



SS+ = 1 



(31) 



In our three terminal geometry nonlocal conductance 
arises only from trajectories that cross both interfaces, 
as illustrated in Fig.[4j Accordingly, the above boundary 
conditions should be employed at both NS interfaces. 

Finally, one needs to specify the asymptotic bound- 
ary conditions far from NS interfaces. Deep in metallic 



electrodes we have 



„K 



7f 


= 7f 


= lt = 


it - o, 




(32) 


h (s + 


eVi), 


if = 


-h (s - 


eV x ), 


(33) 




= 7 2 * 


-it- 


it = o, 




(34) 


ho(e + 


eV 2 ), 


x 2 = 


-h (e - 


eV 2 ), 


(35) 



where ho(e) — tanh(e/2T) - equilibrium distribution 
function. In the bulk of superconducting electrode we 
have 



jy = -a(e)ia 2 , ly = a*(e)icr 2 , (36) 

if = -[l-\a(e)\ 2 ]h (e), (37) 

j§ = a(e)icr 2 , ^,=-a*{e)ia 2l (38) 

= {l-\a(e)\ 2 }h (e), (39) 



where we denoted a(e) — — (e — \J e 2 — A 2 )/A. 



D. Green functions 

With the aid of the above equations and boundary con- 
ditions it is straightforward to evaluate the quasiclassi- 
cal Green-Keldysh functions for our three-terminal device 
along any trajectory of interest. For instance, from the 
boundary conditions at the second interface we find 



T§ = ia(e)A 2 a 2 , 



(40) 



where A 2 = S 2 ' 2 '<t 2 S_ 2 , 2 ,<t 2 . Integrating Eq. (JTJ) along the 
trajectory connecting both interfaces and using Eq. (j4"0|) 
as the initial condition we immediately evaluate the Ric- 
cati amplitude at the first interface: 



R aA 2 + (aA 2 e + A)Q 
lv=l l-(aA 2 A + s)Q °* 
tanh [iflL/vF 



(41) 



n = \Js 2 - A 2 . (42) 



Employing the boundary conditions again we obtain 

rf = iSivK 21 - [aA 2 + (aA 2 e + A)Q] cr 2 S+,, (43) 

ff = -iaS+, x a 2 SvvK^ [1 - (aA 2 A + e)Q] S v \,S V i, 

(44) 

where 

Kij = (1 - a 2 AiAj) - [e(l + a 2 AiAj) + Aa(A t + Ajj\ Q, 

(45) 

A 1 =a 2 S+ v <j 2 S vv . (46) 

We also note that the relation {Y R ^ A )+ = f A < R and 
(yR>^)+ = ZjA,R j^ajjgg ft unnecessary (while redundant) 
to separately calculate the advanced Riccati amplitudes. 



5 




FIG. 4: Riccati amplitudes for incoming and outgoing trajectories for an NSN structure with two barriers. The arrows define 
quasiparticle momentum directions. We also indicate relevant Riccati amplitudes and distribution functions parameterizing 
the Green-Keldysh function for the corresponding trajectory. 



Let us now evaluate the distribution functions at both 
interfaces. With the aid of the boundary conditions at 
the second interface we obtain 



X* = S 2 > 2 'S+ 2 , (1 - \a\ 2 ) h (e) + S 2 , 2 S+ 



r K - 

2 Jj 2 



- \a\'S 2>2 ,a 2 St 2 ,S 22 ,a 2 S+ 2 ,x«. (47) 

Integrating Eq. © along the trajectory connecting both 
interfaces with initial condition for we arrive at the 
expression for Xy 



-i vK, 



xf, = [1- (aA 3 A + e)Q]~ 1 X. 

x (1 - tanh 2 iLn/v F ) [1 - (aA 2 A + e)Q] 



(48) 



Then we can find distribution functions at the first in- 
terface. On the normal metal side of the interface we 
find 

X? = rgsfrjf + t? lX «t« + + a?^ (l - l«| 2 ) h (s) 

(49) 

where 

r£ = Su.Kn 1 [(1 - (aA 2 A + e)Q)S+ v S+~ 
-a(aA 2 + (aA 2 e + A)Q)a 2 Sf, 1 ,a 2 S^ 1 

tg = S 1V K£{1 - (aA 2 A + e)Q), 
af z = iS iy K 21 1 (aA 2 + (aA 2 e + A)Q)a 2 Sf, v 



(50) 

(51) 
(52) 



The corresponding expression for X^ is obtained analo- 
gously. We get 

X? = fOsffS* - (1 - |a| 2 ) h (e) - a%x$ag . 

(53) 



where 



S v \S vv n 2 {\ - (aA 2 A + e)Q)- 



(54) 



-S-,\a 2 S vv a(aA 2 + (aA 2 e + A)Q) K^a 2 St v , 

ill = Sf^S+^a^l - (aA 2 A + e)Q)K^a 2 S+ n ,, (55) 
an = iaS+^SyyK^il - (aA 2 A + e)Q). (56) 

Combining the above results for the Riccati amplitudes 
and the distribution functions we can easily evaluate the 
Keldysh Green function at the first interface. For in- 
stance, for the trajectory l out (see Fig. @| we obtain 



9 U 



(57) 



The Keldysh Green function for the trajectory l m is eval- 
uated analogously, and we get 



g? in = 2x«, g$ n =2(X? -rfxfrf 



(58) 



III. NONLOCAL CONDUCTANCE 
A. General results 

Now we are ready to evaluate the current Ii across the 
first interface. This current takes the form: 



Go 
8e 



h = I? TK {Vi) I de Sp(f 8 $£.« - t 3 g%), (59) 



where 



Gq — 



87172A/1W2 



(60) 



is the normal state nonlocal conductance of our device 
at fully transparent interfaces, Pf7i(2) is normal to the 
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first (second) interface component of the Fermi momen- 
tum for electrons propagating straight between the inter- 
faces, N\.2 = PpAxp/Aif define the number of conduct- 
ing channels of the corresponding interface, R q = 27r/e 2 
is the quantum resistance unit. 

Here lf TK (Vi) stands for the contribution to the cur- 
rent through the first interface coming from trajectories 
that never cross the second interface. This is just the 
standard BTK contribution 2 ^. The non-trivial contri- 
bution is represented by the last term in Eq. (|59[) which 
accounts for the presence of the second NS interface. We 
observe that this non-local contribution to the current is 
small as oc l/p F L 2 (rather than oc l/p F L 3 as suggested 
in Ref. This term will be analyzed in details below. 

The functions gf in and gfLt are the Keldysh Green 
functions evaluated on the trajectories l m and l out re- 
spectively. Using the above expression for the Riccati 
amplitudes and the distribution functions we find 

Sp(T 3 <^Lt -r 3 <^L) = 

= 2 Sp[r>* + - f f f f + - l](h (s + eVj) - h (e))- 
- 2S P [f3fg + - Tfrf + - l](h (e - eVx) - h (s)) + 
+ 2(1 - tanh 2 iLQ/v F )x 
x Sp[K^{S V2 S+ 2 (h (e + eV 2 ) - h Q (e))+ 
+ \a\ 2 S 2 '2'a 2 S+ 2 ,S 22 , ( 7 2 S+ 2 ,{ho(e-eV 2 )-ho(e))}K+; 1 x 
x (S^Sw - W^S+^S^.S+^Svv)}, (61) 

where we explicitly used the fact that in equilibrium 
Sp(f 3 <?fL - f 3 <?£„) = 0. Substituting (JSTJ) into 1(55)1. 
we finally obtain 

h = Ii TK (Vi) + hi(Vi) + I 12 {V 2 ). (62) 

The correction to the local BTK current (arising from 
trajectories crossing also the second NS interface) has 
the following form 

S P [?M + - fff f - l](h (s + eVx) - h (e))- 

- Sp[f*f« + - Tfrf + - l](h (s - eVx) - ho(e))}, 

(63) 

while for the cross-current we obtain 

Iu(Vn) = —^ J de{\- t&nh 2 im/v F )x 

x Sv[K^{S 2 , 2 S+ 2 {h (e + eV 2 ) - h Q {e))+ 
+ \a\ 2 S 2 , 2 ,a 2 S+ 2 ,S 22 ,n 2 S+ 2 ,(ho{e-eV 2 )-h a {e))}K+; 1 x 
x (S+.Su, - W^S+^S.^S+^Srr)}. (64) 

Eqs. ((62)) - (l64)) fully determine the current across the first 
interface at arbitrary voltages, temperature and spin- 
dependent interface transmissions. 



i(T2a(e) 




FIG. 5: Diagrams representing four different contributions 
to the cross-current lis ((M)) . Solid (dotted) lines correspond 
to propagating electron-like (hole-like) excitations and <2i = 
/ cosh(iLQ/v F ). 



In right hand side of Eq. fM)) we can distinguish four 
contributions with different products of S'-matrices. Each 
of these terms corresponds to a certain sequence of ele- 
mentary events, such as transmission, reflection, Andreev 
reflection and propagation between interfaces. Diagram- 
matic representation of these four terms is offered in Fig. 
(5) The amplitude of each of the processes is given by the 
product of the amplitudes of the corresponding elemen- 
tary events. For instance, the amplitude of the process 
in Fig. 5c is / = —iSii>t 2 iS 2 > 2 >aa 2 S_ 22 ,- In Eq. (|E4f this 
process is identified by the term Sp(// + ) with the hole 
distribution function as a prefactor. It is straightforward 
to observe that the processes of Fig. 5a, 5b and 5d corre- 
spond to the other three terms in ((64)) . We also note that 
the processes of Fig. 5a and 5d describe direct electron 
(hole) transport, while the processes of Fig. 5b and 5c 
correspond to the contribution of CAR. 

Assuming that both interfaces possess inversion sym- 
metry as well as reflection symmetry in the plane normal 
to the corresponding interface we can choose 5-matrices 
in the following form 



— Sri' — Sji — Sjn> 



1%) 







3 -i«i/2 



U+(<p), (65) 



£%)t 




- tfl /2 ) u (66) 



and 



S 22 — S 2 > 2 > — S_ 22 — S_ 2 , 2 , — 

l~R^e l ^l 2 



5*2'2 




-i8 2 /2 



, (67) 



(68) 
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Here Ri(2)Ul) = 1 — ^i(2)T(l) are the s P m dependent 
reflection coefficients of both NS interfaces, O12 are spin- 
mixing angles and U (<p) is the rotation matrix in the 
spin space which depends on the angle <p between polar- 
izations of the two interfaces, 

rrf \ 1 ■ /i\ ( cos(W2) -ism(ip/2)\ 

(69) 

Before turning to a detailed calculation of the electric 
current let us briefly address the issue of the spin current 
conservation. It is worth pointing out that in general the 
spin current needs not to be conserved in heterostruc- 
tures with spin active interfaces, see, e.g., Ref. and 
further references therein. Only in certain specific sit- 
uations such conservation can take place. For instance, 
one can easily check that a single barrier with S'-matrix 
(f6"T)) -(f6"5 |) preserves the spin current conservation^. Mak- 
ing use of the general expressions for the Green-Keldysh 
functions we have verified that in our two barrier struc- 
ture with interface ^-matrices (|55|) - (|55| the spin current 
is in general not conserved even in the normal state. 



For instance, spin accumulation on the first barrier is 
controlled by the combination (D21 — D2i)simp which 
vanishes only for collinear barrier polarizations or spin- 
isotropic second interface. We see no particular rea- 
sons to expect conservation of the spin current in gen- 
eral NSN structures with non-collinear interface polariza- 
tions. However, this general case requires a more detailed 
analysis which goes beyond the frames of the present pa- 
per and will be published elsewhere. 

Now we turn to the analysis of the electric current. 
Substituting the above expressions for the S'-matriccs 
into Eqs. (|6"3")) and we arrive at the final results for 
both iii(Vi) and Ii 2 (V 2 ) which will be specified further 
below. 



B. Cross-current 

First let us consider the cross-current 112(^2). From 
the above analysis we obtain 



de 



tanh — — — - — tanh — 
2T 2T 



1 — tanh 2 iLfl/vp 



W(zi, z 2 , e, ip) 

{ [D H D 2l - \a\ 2 D 1] D n {R li + R 2] ) + | a | 4 J D u i?i T C 2i i?2 T ] \K(z u z 2 , e)\ 2 cos 2 (^/2) + 
+- [Di T D 2T - \a\ 2 D n D 2] (R^+R 2i ) + \a\ 4 D^R n D 2 ^R 2l ] \K{z{, z*, e)\ 2 cos 2 (^/2)+ 
+ [D l1 D 2l - |a| 2 Z? u Z? 2i (i? 1T + i? 2T ) + |a| 4 J Di T i? u I? 2i .R 2T ] \K(z* 1 ,z 2 ,e)\ 2 sin 2 (y>/2)+ 
■f [D U D 2T - \a\ 2 D 1] D 2] (R n +R 2i ) + \a\ 4 D n R n D 2] R 2i ] \K(z u z*, s)\ 2 sm 2 (^/2)}, 



(70) 



where we define 



K(zi,z 2 , s) = (1 - a 2 ziz 2 ) — [e(l + a 2 ziz 2 ) + Aa(z 1 + z 2 )] Q, 
W( Zl ,z 2 , e, if) = \K(z u Z2, e)K{zl,z*,e) cos 2 (^/2) + K(zt,z 2 , e)K{ Zl ,z*, e) sin 2 (^/2)| 2 



(71) 
(72) 



and Zi = ^jR^Rii exp(i6i) ( i — 1,2). 



r 



Eq. ((70|) represents our central result. It fully deter- 
mines the non-local spin-dependent current in our three- 
terminal ballistic NSN structure at arbitrary voltages, 
temperature, interface transmissions and polarizations. 

Let us introduce the non-local differential conductance 



G 12 (V 2 



dh 
dV 2 



dh 2 (V 2 ) 
dV 2 • 



(73) 



Before specifying this quantity further it is important 
to observe that in general the conductance Gi 2 (V 2 ) is 
not an even function of the applied voltage V 2 - This 
asymmetry arises due to formation of Andreev bound 
states in the vicinity of a spin-active interfac e) 35 ' 36 . It 



disappears provided the spin mixing angles 9\ and 9 2 
remain equal to or tt. 

In the normal state we have I\ 2 {V 2 ) = — Gn 12 V2, where 

Gn 12 = ^ [(D H D 2l + Di T D aT ) cos 2 (^/2) + 

+ (D 1 ^D n +D n D 2 ^) sin 2 (ip/2)]. (74) 

Turning to the superconducting state, let us consider 
the limit of low temperatures and voltage T, V2 <C A. 
In this limit only subgap quasiparticles contribute to the 
cross-current and the differential conductance becomes 
voltage-independent, i.e. I± 2 — —Gi 2 V 2 , where 
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G12 = G (l - tanh 2 LA/v F )< 



\K( Zl ,z 2 , 0)| 2 cos 2 ((^/2) + \K( Zl , z*, 0)| 2 sin 2 (^/2) 



m n w n n , |g(z 1 ,z 2 ,0)| 2 co S 2 (^/2)-|X(z 1 ,z 2 * ; 0)| 2 sin 2 (y/2) 1 

(L>n - L>ii )(L>2T _ ^21) o >■ (75) 

z 2 , 0)P cos^/2) + \K( Zl , z*, 0)P sin 2 (^/2)) 2 J 



0.25 




At the first sight, this result might appear counterintu- 
itive since the behavior of ordinary (local) AR is just the 
opposite: It reaches its maximum at full barrier transmis- 
sions. The physics behind vanishing of CAR for perfectly 
trasparent NS interfaces is simple. One observes (cf. Fig. 
1) that CAR inevitably implies the flow of Cooper pairs 
out of the contact area into the superconducting termi- 
nal. This flow is described by electron trajectories which 
end deep in the superconductor. On the other hand, 
it is obvious that CAR requires "mixing" of these tra- 
jectories with those going straight between two normal 
terminals. Provided there exists no normal electron re- 
flection at both NS interfaces such mixing does not occur, 
CAR vanishes and the only remaining contribution to the 
non-local conductance is one from direct electron transfer 
between N-terminals. 



eV 2 /A 



FIG. 6: Zero temperature differential non-local conductance 
as a function of voltage at zero spin-mixing angles #1,2 = 0. 



In the case of spin-isotropic interfaces Eqs. (|75|) and (|70|) 
reduce to our previous results^ 2 -. 

Provided at least one of the interfaces is spin-isotropic, 
the conductance ((75)) is proportional to the product of 
all four transmissions D^D^Z^t-^Xj i-e- it differs from 
zero only due to processes involving scattering with both 
spin projections at both NS interfaces. As in the case 
of spin-isotropic interfaces^ the value Gi 2 (|75|) gets 
strongly suppressed with increasing L, and at sufficiently 
high interface transmissions this dependence is in general 
non-exponential in L. In the spin-degenerate case for a 
given L the non-local conductance reaches its maximum 
for reflectionless barriers D12 — 1. In this case we arrive 
at a simple formula 

Gvi = G (l - tanh 2 LA/v F ). (76) 

We observe that for small L <C vp/A the conductance 
G12 identically coincides with its normal state value 
Gjv 12 = Go at any temperature and voltageA 2 ". This re- 
sult implies that CAR vanishes for fully open barriers. 
Actually this conclusion is general and applies not only 
for small but for any value of L, i.e. the result (|76")l is de- 
termined solely by the process of direct electron transfer 
between N-terminals for all L. 



This situation is illustrated by the diagrams in Fig. 
[5] It is obvious that in the case of non-reflecting NS in- 
terfaces only the process of Fig. 5a survives, whereas all 
other processes (Fig. 5b, 5c and 5d) vanish for reflection- 
less barriers with i? 1 ( 2 )|(|) = 0. The situation changes 
provided at least one of the transmissions is smaller than 
one. In this case scattering at SN interfaces mixes up tra- 
jectories connecting Ni and N2 terminals with ones going 
deep into and coming from the superconductor. As a re- 
sult, all four processes depicted in Fig. 5 contribute to 
the cross-current and CAR contribution to G12 does not 
vanish. 



Let us also note that the statement about the ab- 
sence of CAR for highly transparent interfaces was in- 
dependently made in a recent paper—. Although no 
derivation supporting this statement was presented, this 
statement^! appears to be based on the BTK-like de- 
scription of strictly Id NSN structures. According to our 
general analysis, no processes presented in Figs. 5b, 5c 
and 5d would be possible in that case and, hence, CAR 
should be prohibited in Id NSN systems for any transmis- 
sion. This observation, however, can be easily overlooked 
if one only deals with the solutions of the Bogolyubov- 
de Gennes equations (equivalent to finding retarded and 
advanced Green functions) and does not directly eval- 
uate the electron distribution function (or the Keldysh 
function) in the S-terminal. 
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0.4 



0.3 



0.2 



0.1 



-0.1 



tp = 
tp — 7r/4 
<p = tt/2 
tp — 3tt/4 

tp = 7T 



\ 

\ \ 



/ £>i T = 0.2 
£>u = 0-5 



L = vf/A 
.(, 2 \. 0! = vr/2 
= vr/4 

-0.3 



V..-7 



D 2T = 0.3 - 
D 2I = 0.7 



-2 -1.5 I -0.5 0.5 
eV 2 /A 



1.5 



FIG. 7: The same as in Fig. 6 for 9i = vr/2, 2 = vr/4. 



In the limit \eV2\1T <C A and at zero spin- mixing an- 
gles 9\ ,2 =0 from Eq. 1(75]) we obtain 

_ 1 - tanh 2 LA/v F , 
Gl2 ~ G ° 1^,^2,0)12 1^T^U^2T^2| + 

+ ( J D 1T - J D u )( J D 2T - J D 2i )cos^}. (77) 

In the lowest (first order) order in the transmissions of 
both interfaces and for collinear interface polarizations 
Eq. (|77|) reduces to the result by Falci et aL— pro- 
vided we identify the tunneling density of states TVo-Dif, 
N Dii, iVo-C^ti an d N D 2 i with the corresponding spin- 
resolved densities of states in the ferromagnetic elec- 
trodes. For zero spin-mixing angles and low voltages the 
L-dependence of the nonlocal conductance Gi 2 reduces 
to the exponential form G12 oc exp(— 2LA/vf) either in 
the limit of small transmissions or large L ^> vf/A. 

At arbitrary voltages and temperatures the cross- 
current has a simple tp dependence in the limit of zero 
spin mixing angles (61,2 = 0) 



/l2(^2)=/l 2 (^ = 0,y 2 )cOS 2 



(p/2)- 



/ 12 (^ = 7r,y 2 )sin 2 ((^/2), (78) 



i.e. in this limit at any tp the nonlocal current is equal 
to a proper superposition of the two contributions cor- 
responding to parallel (tp = 0) and antiparallel (tp = it) 
interface polarizations. Some typical curves for the dif- 
ferential non-local conductance are presented in Fig. [6] 
at sufficiently high interface transmissions and zero spin 
mixing angles 612 = 0. 

Let us now turn to the limit of highly polarized inter- 
faces which is accounted for by taking the limit of vanish- 
ing spin-up (or spin-down) transmission of each interface. 
In this limit our model describes an HSH structure, where 
H stands for fully spin-polarized half-metallic electrodes. 



In this case we obtain (-Dij 
and D21 = 0) 



£>i, D n =0,D 



21 



Do 



IM) 



Go 



de [ho(e + eV 2 ) - h Q (e)] x 



1 — tanh 2 iLVI/vf 

x — ttt-, — — -Dl-D 2 X 



W(zi,z 2 ,e,tp) 
x {[l + |a| 4 ] \K(zlz* 2 ,s)\ 2 C os 2 (tp/2)- 

-2|a| 2 |X(z 1 ,z 2 *, £ )| 2 sin 2 (^/2)}. (79) 

We observe that the nonlocal conductance has opposite 
signs for parallel (tp — 0) and antiparallel (tp = it) in- 
terface polarizations. We also emphasize that, as it is 
also clear from Eq. (|77|) . the cross-conductance Gi 2 of 
HSH structures - in contrast to that for NSN structures 
- does not vanish already in the lowest order in barrier 
transmissions Dij_D 2 j. 

In general the non-local conductance is very sensitive 
to particular values of the spin-mixing angles 9\ and 62, 
as illustrated, e.g., in Fig. [7] Comparing the voltage 
dependencies of the nonlocal conductance evaluated for 
the same transmissions and presented in Figs. [S] and [7J 
we observe that they can differ drastically at zero and 
non-zero values of #i i2 . 

At low voltages and temperatures and at zero spin mix- 
ing angles the non-local conductance of HSH structures 
is determined by Eq. (J77J with Di± = D 2 i = 0. For fully 
open barriers (for "spin-up" electrons) Dif = D 2 -[ = 1 
we obtain 



Gi 2 = Go(l — tanh 2 LA/vp) cos tp. 



(80) 



Interestingly, for tp — this expression exactly coincides 
with that for fully open NSN structures, Eq. (76|). At 
the same time for small L the result (|80|) turns out to 
be 2 times bigger that the analogous expression in the 
normal case, i.e. for fully open HNH structures, cf. Eq. 
([71)1 . This result can easily be interpreted in terms of 
diagrams in Fig. 5. We observe that - exactly as for 
the spin degenerate case - CAR diagrams of Fig. 5b, c 
vanish for reflectionless barriers, whereas diagrams of Fig. 
5a, d describing direct electron transfer survive and both 
contribute to Gi 2 . Thus, CAR vanishes identically also 
for fully open HSH structures. The factor of 2 difference 
with the normal case is due to the fact that the diagram 
of Fig. 5d vanishes in the normal limit. 

Finally, let us compare our general results, Eq. ([70]) 
and below, with the the analogous results^. In order 
to account for higher order tunneling events in the case 
of symmetric interfaces the authors^ employed effec- 
tive "dressing" of both EC and CAR contributions by 
all higher order local Andreev processes at both inter- 
faces. At the same time, all higher order non-local pro- 
cesses were ignored in Ref. ITol . Our analysis - which in- 
cludes all processes to all orders - demonstrates that this 
approximation^ may be appropriate only in the limit 
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L 3> vp/A in which case non-local effects are exponen- 
tially suppressed. Otherwise higher order non-local pro- 
cesses remain important and need to be fully accounted 
for. 

Comparing our exact results and those of Ref. [13 (with 
omitted higher order non-local processes) we observe that 
the disagreement between them grows with increasing 
transmissions and becomes maximal at full transmissions 
and small values of L. The most significant differences 
are: (i) the dependence of G\ 2 on L is non-exponential 
whereas in Ref. it remains exponential at all L, (ii) at 
small L and high transmissions the linear conductance 
G12 evaluated in Ref. 03 (see e.g. Eqs. (33), (34) of 
that work) turns out to be several times smaller than our 
results in both cases of NSN and HSH structures and (iii) 
taking into account only local Andreev reflection events 
is not sufficient to correctly account for vanishing CAR 
in the limit of fully open interfaces. 

C. Correction to BTK 

Using the above formalism one can easily generalize 
the BTK result to the case of spin-polarized interfaces^!. 



For the first interface we have 



lBTK {Vi) = p_ J d£[ho{£ + eVl )- ho{£)]{ l + \ a \2) x 



Here transmission and reflection coefficients as well as the 
spin mixing angle depend on the direction of the Fermi 
momentum. In the spin-degenerate case the above ex- 
pression reduces to the standard BTK result 2 -. 

Evaluating the nonlocal correction to the BTK current 
due to the presence of the second interface we arrive at 
a somewhat lengthy general expression 



hx{Vx) = ^- / de{h (e + eV,) - h Q (e))——^ -{ 2W( Zl , z 2 , e,tp)- 

2e J W(zi,Z2,e,(p) I 

~R 11 \cos^^/2)K{z 1 /R lh z 27 s)K(zlz; 7 e)+sin 2 { V /2)K(z 1 /R lh z;,e)K{z* l7 z 2 ,e)\ 2 - 
-R n \cos 2 (v/2)K{zl/R n ,z* 2 , e)K(zi,z 2 , e)+ sin 2 ( l p/2)K(z* 1 /R ll ,z 2 , e)K( Zl ,z*, e)\ 2 } + 

+ ^ [ de(h (e + e Vl ) - h (e)) { 
4e / Wlz-i . zo. e, w) I 



4e J ' W(zi,z 2 ,s,tp) 

\a\ 2 \cos 2 ((p/2)K(Q, z 2 , s)K{zl,z 2 , e)+ sin 2 (^/2)if(0, z 2 , s)K(z*, z 2 , e) | 2 + 
+ |a| 2 |cos 2 (^/2)K(0, z*, e)K( Zl , z 2 , e)+ sin 2 (<p/2)K(0, z 2 , e)K{z u z* 2: e) | 2 + 

+ H2 |cos 2 ^/2)^'(z 2 * ,e)K(z u z 2 , e)+ sm 2 (<p/2)K'(z 2 ,e)K(z 1 , z*,e)\ 2 + 
\a\ 

+ -\r\cos 2 (ip/ 2) K'(z 2 , e)K(zl,z* 2 , e)+ sin 2 {ip/2)K '{z*, e)K {z{, z 2 , e)\ 2 } + 
M " J 

+ ^ J R 2T i? 2i sin 2 (0 2 /2)sin 2 (^/2)cos 2 ^/2) / de(h (e + eVi) - h (e)) {1 ~ ^ ' 



W(zi,z 2 ,e,ip) 

x [|a| 2 (i? 2 T + D 2 n ) - 2\a\ i D 1] D n {R 1] + R n ) + \a\ e (D^R 2 n + D 2 ^)} , (82) 

I 

where K'(z 2 ,e) — dK(z\, z 2 , e)/dz\. This expression mixing angles 6*1,2 = in which case we obtain 
gets significantly simplified in the limit of zero spin- 



hi{Vi) = J de{h (e + eV 1 



\K{ Zl /R lh z 2 ,e)\ 2 \K( Zl /R n ,z 2 ,e) 
\K(z u z 2 ,e)\* 1 \K{z u z 2 ,e)\* 



if) r) He)\ 2 \K(0,z 2 ,e)\ 2 + \K>(z 2 ,e)\ 2 /\a(e)\ 2 \ 
1T U \K( Zl ,z 2 ,eW j' 

(83) 
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In contrast to the expression for the cross-current I\ 2 
(cf. Eq. (|78f ). in the limit of zero spin-mixing angles 
the correction I\\ to the BTK current does not depend 
on the angle ip between the interface polarizations. In 
particular, at |eVi|,T <C A we have In = G\\V\ where 



Gn — Go(Dif+Dii 



(1 - zf)(l - tanh 2 LA/vp) 
1 + ziz 2 + {z\ + z 2 ) tanhLA/^i?] 
(1 + z 2 tanhiA/w F ) 2 + 3(z 2 + tanhLA/w F ) 



[1 + z\z 2 + (zi + z 2 ) tanhLA/w^] 



In the tunneling limit -Dif, Dij, D 2 ^, D 2 ± < 1 we repro- 
duce the result of Ref. d 

Gn = ^-(D 1 - t + D n )(D v +D 2l )e X p(-2LA/v F ), (85) 

which turns out to hold at any value (p. 

As compared to the BTK conductance the CAR cor- 
rection (f8"2")l contains an extra small factor A 2 /L 2 and, 
hence, in many cases remains small and can be neglected. 
On the other hand, since CAR involves tunneling of one 
electron through each interface, for strongly asymmetric 
structures with D\^,Dxi <C 1 and D 2 ^,D 2 i ~ 1 it can 
actually strongly exceed the BTK conductance. Indeed, 
for Dm <C 1, R 2 -[R 2 i <C 1 and provided the spin mixing 
angle B\ is not very close to 7r from Eq. (|82|) we get 



Gn — 



G (D 



n 



D 



n) 



i.e. for 



the contribution 



cosh(2LA/u F ) + cos 6»i smh(2LA/vF) 



^^ y <^exp(-2LA/, F ) 



(86) 



three-terminal structures with spin- active interfaces. 
Our theory is based on the quasiclassical formalism 
of energy-integrated Green-Eilenberger functions supple- 
mented by appropriate boundary conditions describing 
spin-dependent scattering at NS interfaces. Our ap- 
proach applies at arbitrary interface transmissions and 
allows to fully describe non-trivial interplay between 
^^spin-sensitive normal scattering, local and non-local An- 
2 dreev reflection at SN interfaces. Our main results are 
the general expressions for the non-local cross-current 
Ii2, Eq. ([70| . and for the non-local correction In to 
(84ihe BTK current, Eq. (JS2J). These expressions provide 
complete description of the conductance matrix of our 
three-terminal NSN device at arbitrary voltages, tem- 
perature, spin-dependent transmissions of NS interfaces 
and their polarizations. 

Our analysis allows to predict and analyze a rich va- 
riety of interesting properties of such structures. One of 
our predictions is that in the case of ballistic electrodes 
no crossed Andreev reflection can occur in both NSN 
and HSH structures with fully open interfaces. Beyond 
the tunneling limit the dependence of the non-local con- 
ductance on the size of the S-electrode L is in general 
non-exponential and reduces to G12 oc exp(— 2LA/vp) 
only in the limit of large L. For hybrid structures half- 
metal-superconductor-half-metal we predict that the low 
energy non-local conductance does not vanish already in 
the lowest order in barrier transmissions G12 oc D\^D 2 ^. 
These and other our predictions can be directly tested in 
future experiments on NSN hybrid structures, including 
systems with ferromagnetic and half-metallic electrodes. 



(D 



BTK 



may well exceed the BTK term 
Di^Dii- The existence of such a non-trivial 



regime further emphasizes the importance of the mech- 
anism of non-local Andreev reflection in multi-terminal 
hybrid NSN structures. 



IV. CONCLUSIONS 

In this paper we developed a non-perturbative the- 
ory of non-local electron transport in ballistic NSN 
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